Astrometry in the Age of the Next Generation of Large Telescopes 

ASP Conference Series, Vol. 000, 2005 

P. Kenneth Seidelmann and Alice K. B. Monet 

Science Results Enabled by SDSS Astrometric 
Observations 

Zeljko Ivezic 1 , Nicholas Bond 2 , Mario Juric 2 , Jeffrey A. Munn 3 , Robert 
H. Lupton 2 , Jeffrey R. Pier 3 , Gregory S. Hennessy 4 , Gillian R. Knapp 2 , 
James E. Gunn 2 , Constance M. Rockosi 5 , Tom Quinn 1 

1 University of Washington, 2 Princeton University, ^USNO Flagstaff, 
4 USNO Washington D.C., ^University of California-Santa Cruz 

Abstract. We discuss several results made possible by accurate SDSS astro- 
metric measurements in a large sky area, with emphasis on asteroids and stellar 
proper motions obtained by comparing POSS and SDSS. SDSS has observed 
over 200,000 moving objects in five photometric bands, corresponding to about 
two orders of magnitude increase over previous multi-color surveys. These data 
were used to extend the measurement of asteroid size distribution to a smaller 
size limit, to demonstrate that asteroid dynamical families, defined as clusters 
in orbital parameter space, also strongly segregate in color space, and to dis- 
cover a correlation between asteroid age and colors. A preliminary analysis of 
SDSS-POSS proper motions for ~1 million M dwarf stars demonstrates that, in 
the 0.1-1 kpc distance range, the rotational velocity and its dispersion for disk 
stars increase with the distance from the Galactic plane. 



1. Introduction 

The astrometric obervations obtained by modern massive digital sky surveys, 
such as SDSS, 2MASS, and FIRST, are enabling studies that have not been 
possible before. In this contribution we describe several science results enabled 
by SDSS astrometric observations, with emphasis on observations of asteroids 
and proper motions obtained by comparing POSS and SDSS astrometric mea- 
surements. 

The Sloan Digital Sky Survey (SDSS) is currently mapping one quarter of 
the sky in five optical bands (ugriz, Fukugita et al. 1996; Gunn et al. 1998) to a 
depth of r ~ 22.5, accurate to 0.02 magnitudes (both absolute calibration, and 
root-mean-square scatter for sources not limited by photon statistics; Ivezic et 
al. 2004). Astrometric positions are accurate to better than 0.1" per coordinate 
(rms) for sources with r < 20.5 (Pier et al. 2003), and the morphological informa- 
tion from the images allows reliable star /galaxy separation to r ~ 21.5 (Lupton 
et al. 2002). The survey's coverage of ~ 10 4 deg 2 in the North Galactic Cap and 
of ~ 200 deg 2 in the Southern Galactic Hemisphere will result in photometric 
and astrometric measurements for over 10 8 stars and a similar number of galax- 
ies. Additionally, SDSS will obtain spectra for over 10 6 objects, including 10 6 
galaxies and 10 5 quasars. The recent third public Data Release (DR3) includes 
imaging data for 5282 deg 2 of sky, and catalogs 1.4 x 10 8 objects (Abazajian et 
al. 2005). A detailed technical description of SDSS is given by Stoughton et al. 
(2002). 
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The SDSS astrometric reductions, and the measurements of their success, 
are described in detail by Pier et al. (2003). Briefly, the SDSS absolute astro- 
metric accuracy is better than 100 mas, with relative (band-to-band) accuracy 
of about 20-30 mas (rms, for sources not limited by photon statistics). In ad- 
dition to providing positions for a large number of objects with a remarkable 
accuracy (and thus enabling recalibration of other less accurate surveys, as de- 
scribed below), an important characteristic of SDSS astrometric observations is 
that measurements in five photometric bands are obtained over a five minute 
long period (with 54 sec per exposure). The multi-color nature allows the dis- 
covery of the so-called Color Induced Displacement binary stars (for details see 
Pourbaix et al. 2004), and the time delay allows the recognition of moving 
objects, discussed next. 

2. SDSS Observations of Solar System Objects 

Although primarily designed for observations of extragalactic objects, the SDSS 
is significantly contributing to studies of solar system objects, because asteroids 
in the imaging survey must be explicitly recognized to avoid contamination of 
the quasar samples selected for spectroscopic observations (Lupton et al. 2002). 
The SDSS has already increased the number of asteroids with accurate five-color 
photometry by about two orders of magnitude (to over 200,000), and to a limit 
more than five magnitudes fainter than previous multi-color surveys (Ivezic et 
al. 2001). 

About 43,000 of those 200,000 objects have been associated with previously 
known asteroids that have well determined orbital elements (Juric et al. 2002). 
Both SDSS data and orbital elements are available from the public SDSS Mov- 
ing Object Catalog^ (Ivezic et al. 2002a). The SDSS observations of objects 
with known orbits vividly demonstrate that asteroid dynamical families, defined 
as clusters in orbital parameter space, also strongly segregate in color space 
(Ivezic et al. 2002b). This segregation indicates that the variations in chemical 
composition within a family are much smaller than the compositional differ- 
ences between families, and strongly support earlier suggestions that asteroids 
belonging to a particular family have a common origin. 

Asteroid colors measured by SDSS are well correlated with the family age 
(Jedicke et al. 2004), which provides a direct evidence for space weathering, and 
offers a method to date asteroids using their SDSS colors. The effects of space 
weathering can also explain the color variability of asteroids measured by SDSS 
(Szabo et al. 2004). 

3. Proper Motions Determined from SDSS and POSS Observations 

Munn et al. (2004, hereafter M04) have recently presented an improved proper- 
motion catalog for ~8 million stars, based on combining the USNO-B and SDSS 
catalogs in the 2099 deg 2 large area of sky covered by SDSS Data Release 1. 
USNO-B positions are recalibrated using SDSS galaxies, which results in smaller 



1 Available from http: / /www.sdss.org/dr2 /products/value_added / index.html 
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Figure 1. An illustration of the differences in the distribution of point 
sources in SDSS color-magnitude and color-color diagrams induced by re- 
quiring a large proper motion. The distribution of all point sources is shown 
by contours, and two subsets of sources with large proper motion are shown 
by symbols (small: 4-7 mas/yr; large: 7-15 mas/yr). The white dwarfs are 
easily recognized by their blue u — g color (u — g < 0.8), while the majority of 
other sources with large proper motions are main sequence stars. M dwarfs, 
at distances up to ~1 kpc, dominate the large proper motion sample. 



astrometric errors and the placement of proper motions on an absolute reference 
frame. Tests based on several tens of thousand quasars spectroscopically con- 
firmed by SDSS show that the proper motion errors are ~3 mas/yr for bright 
(g < 18) sources (rms per coordinate), with substantially smaller systematic 
errors (<0.5 mas/yr, and typically ~0.1 mas/yr), and excellent stability across 
the sky. For fainter sources, the proper motion errors increase to ~4 mas/yr at 
g ~ 19, ~5.7 mas/yr at g ~ 20, and to ~7.5 mas/yr at g ~ 20.5, the faint limit 
for the sample discussed here. For stars at 500 pc, these proper motion errors 
correspond to systematic velocity error of <1 km/s, and to random velocity er- 
rors of 7 km/s at the bright end (without including errors in distance), and 20 
km/s at the faint end. 

This accurate proper motion database offers numerous possibilities for study- 
ing the stellar kinematics, but the analysis of such a large and highly-dimensional 
data set is not trivial. At the very least, the proper motion components are func- 
tions of at least four observables: apparent magnitude, distance (or absolute 
magnitude, measured by e.g. r — i color), and the position on the sky (Fig. [[]). 
Additional stellar parameters, such as metallicity (measured from spectra, or 
estimated from the u — g color), are also expected to play a role. Furthermore, 
the probed distance range is large (out to ~15 kpc with main sequence stars, 
Juric et al. 2005) and thus the dependence of the stellar velocity distribution on 
position in the Galaxy could also be important, as well as the deviations from a 
Gaussian distribution. 

In this contribution we limit our analysis to M dwarfs from the M04 proper 
motion catalog because 1) they dominate the high proper motion sample, and 
2) their distances can be reliably determined from a photometric parallax rela- 
tion. A more comprehensive analysis of that catalog, that also includes three- 
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Figure 2. The bottom left panel shows the distance distribution for a sample 
of 881,913 M dwarfs with SDSS-POSS proper motions. The top panel shows 
the sky distribution for a subsample with distances in the range 100-150 
pc in Aitoff projection of galactic coordinates, and the bottom right panel 
compares the proper motion distribution for this subsample before (solid) 
and after (dashed) correction for solar peculiar motion (see § 13. 1 . [) . 

dimensional velocities for ~100,000 stars with SDSS spectra, will be presented 
by Bond et al. (2005, in prep., hereafter B05). 

We select 881,913 M dwarfs from the M04 catalog by requiring 1.25 < 
g -r < 1.50, 0.6 < r - i < 1.6 and 15 < g < 20.5 (Ivezic et al. 2004), and 
estimate their distances using a photometric parallax relation: 

Mi = 4.0 + 10.86 (r - t) - 10.74 (r - if + 5.99 (r - if - 1.20 (r - if. (1) 

These distances have random errors of ~10%, with comparable systematic errors 
(Juric et al. 2005). The sky and distance distributions of selected stars are shown 
in Fig. [2j The size of this sample and its sky distribution make it well suited for 
a determination of solar peculiar motion, discussed next. 

3.1. Determination of Solar Motion Relative to M dwarfs 

The solar peculiar motion with a magnitude of v©, induces a dipole in the 
proper motion distribution for sources at a distance D, with a magnitude of 
~3.2 (vq/15 km/s) (kpc/D) mas/yr, which is easily measurable with the sample 
discussed here. We follow a procedure outlined by Mihalas & Binney (1981), 
and determine the apex coordinates and the magnitude of implied solar motion 
for subsamples selected in narrow distance bins (Fig. . 

The extrapolation of results shown in Fig. [3] to zero distance gives (1^=68°, 
b J 4=32°, and v^=9 km/s), somewhat different from standard values (1^=51°, 
b J 4=23°, and v^=15 km/s; Mihalas & Binney 1981). However, as evident from 
Fig. [3l all three quantities vary as a function of the subsample distance. We find 
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Figure 3. The galactic coordinates of the solar apex (the top two panels), 
implied proper motion magnitude (bottom left), and corresponding solar ve- 
locity (bottom right), determined as a function of distance using SDSS-POSS 
proper motions for a sample of 881,913 M dwarfs. 



that this effect is caused by the increase of the rotational velocity component 
with the distance from the plane, which breaks the dipole assumption of the 
Mihalas & Binney method. Consequently, the apex is biased toward (1~90°, 
b J 4=0°) direction, with a corresponding increase of v^. An illustration of this 
dependence is discussed next. Additional supporting evidence for this claim, 
which utilizes three-dimensional velocity information, will be described in N05. 

3.2. Kinematics as a Function of Z 

Here we present results for the kinematics of M dwarfs observed towards the 
Galactic anticenter. In this direction, the rotational velocity is dominated by 
the longitudinal velocity component measured from proper motion, and thus 
robust results can be derived without radial velocity information. The top two 
panels in Fig. 2] show the dependence of longitudinal and latitudinal velocity 
components, determined from proper motions, for ~50,000 M dwarfs toward 
1~180 (without including correction for solar peculiar motion, which amounts 
to an overall shift of the longitudinal component by ~10 km/s). While the 
median values of the latitudinal component are assuringly consistent with zero, 
the medians of the longitudinal component increase with the distance from the 
plane, Z, as vj = (12 + 22 |Z/kpc| 1,3 ) km/s (with the uncertainty of ~10% for 
all three fitted parameters). 

From the data presented here it cannot be ruled out that the observed 
variation of v; actually reflects a dependence on the radial distance from the 
Galactic center (because the range of galactic latitudes observed by SDSS to- 
wards 1~180 is limited). However, N05 demonstrate, using three-dimensional 
velocity information, that the Z dependence dominates. Furthermore, they also 
show that the same relation describes the behavior of v; for bluer stars that 
sample distances out to ~3 kpc (though with less reliable distances). 
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Figure 4. The contours in the top left panel show the longitudinal velocity 
component determined from proper motion measurements, vt, as a function 
of distance from the Galactic plane, Z, for 51,978 M dwarfs with 175 < 
I < 185. Toward this direction, the longitudinal velocity is equivalent to 
rotational velocity. The large dots are the median values of vi for a given 
Z bin. The dashed lines represent 2a envelope around the median (a is 
determined using the interquartile range). The line following the medians is 
Vi = (12 + 22 l-Z/kpcl 1 ' 3 ) km/s. The top right panel is analogous, except that 
the latitudinal velocity component is shown. The bottom left panel shows 
the g vs. r — i color-magnitude diagram for a subsample of 9,231 stars with 
250 < Z/pc < 350. The longitudinal velocity component of those stars as 
a function of r — i color is shown in the bottom right panel, with the same 
meaning of symbols and lines as in the top two panels. Note that neither the 
median, nor the distribution width, vary with color, although the apparent 
magnitude, g, varies by ^2 mag. 



3.3. The Shape of the Velocity Distribution 

The analysis described in preceeding sections considers only the first two mo- 
ments of the velocity distribution. Nevertheless, the large number of stars al- 
lows an unprecedentedly robust and accurate analysis of the full shape of the 
velocity distribution. As Figs. [5] and [6] show, the observed shapes significantly 
deviate from Gaussian, and show distortions similar to those expected from the 
asymmetric drift effect (Mihalas & Binney 1981). The shapes also vary with 
distance, a variation that is hard to explain as due to increasing proper motion 
errors toward the faint end (Fig. [6]). The observations of the shape of the veloc- 
ity distribution, and its dependence on the position and stellar tracer, encode 
an enormous amount of information (e.g. Binney & Tremaine 1987) that we are 
only beginning to harvest. 
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Figure 5. The two-dimensional velocity distribution, vi vs. Vb, determined 
from proper motion measurements for two characteristic lines of sight (left: 
175 < I < 185, right: 355 < I < 5) and for the height above the plane in 
the range 100-150 pc. The density of sources is shown by linearly spaced 
isopleths. 



4. Discussion 



The availability of modern massive sky surveys is rapidly changing astronomical 
methodology, and enables studies that were not possible until recently. This is 
particularly true for astrometric data provided by SDSS. A good example is the 
SDSS observations of asteroids: in only a few years, the SDSS has increased the 
sample of asteroids with accurate colors by about two orders of magnitude, and 
enabled detailed and robust studies of correlations between asteroid chemical 
and dynamical properties. 

Another example where SDSS data were crucial is the construction of the 
Munn et al. SDSS-POSS proper motion catalog. This resource is bound to make 
significant contributions to our understanding of the Milky Way kinematics. 
The preliminary results discussed here suggest that these kinematics may be 
significantly more complex than assumed until now. 
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